Vortex coalescence and type-1.5 superconductivity in Sr2Ru04 
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Recently vortex coalescence was reported in superconducting Sr2Ru04 by several experimental 
groups for fields applied along the c-axis. We argue that Sr2Ru04 is a type-f.5 superconductor 
with long-range attractive, short-range repulsive intervortex interaction. The type-1.5 behavior 
stems from an interplay of the two orbital degrees of freedom describing this chiral superconductor 
together with the multiband nature of the superconductivity. These multiple degrees of freedom 
give rise to multiple coherence lengths, some of which are larger and some smaller than the magnetic 
field penetration length, resulting in nonmonotonic intervortex forces. 

PACS numbers: 74.20.De, 74.25.Ha, 74.70.Pq 



The superconducting state of strontium ruthenate is 
of great interest because of experimental evidence that it 
is a chiral spin-triplet superconductor. Support for this 
point of view came through a variety of measurements 
that include: the complete suppression of the supercon- 
ducting transition temperature (Tc) with non-magnetic 
impurities;^ NMR Knight shift measurements that show 
no change in the spin susceptibility with temperature in 
the superconducting phase;^''^ muon spin measurements 
(/iSR) that suggest broken time-reversal symmetry in 
the superconducting state;"* polar Kerr effect showing 
the broken time-reversal symmetry and the presence of 
chirality in the superconducting phase;"'' and phase sen- 
sitive measurements consistent with a chiral spin-triplet 
state. ^'^ While these measurements provide a strong case 
for a chiral spin-triplet superconducting phase, the case 
is not iron clad. In particular, sizable edge currents are 
expected to flow at sample boundaries and at domain 
walls between domains of opposite chirality.® A search 
for these currents has been carried out and they have not 
been observed. ^"^* Additionally, for fields in the basal 
plane, two superconducting phases are predicted^^ and 
the corresponding phase transition between these phases 
has also not been observed. Finally, if spin-orbit coupling 
is sufficiently large, then the spin susceptibility should 
show no change for fields in the basal plane, but is ex- 
pected to decrease with temperature for fields along the 
c-axis. NMR data shows no change on the spin suscep- 
tibility for fields both in the basal plane^ and along the 
c-axis."^ In spite of these puzzles in interpreting the super- 
conducting state as a chiral spin-triplet superconductor, 
the evidence in support of this state remains strong and 
current research is focused on addressing these puzzles 
within this framework. 

Another interesting aspect of the Sr2Ru04 is the multi- 
band nature of its superconductivity. The Fermi surface 
of Sr2Ru04 contains three sheets of cylindrical topol- 
ogy labeled a, /3, and 7.*^ These sheets stem from Ru 
dxz,dyz and Ru dxy orbitals. In particular, the a and (3 
sheets originate from the dxz and dyz orbitals which lead 
to quasi-one-dimensional bands, and the 7 sheet origi- 
nates from the dxy orbitals which leads to a quasi-two- 



dimensional band. It has been shown that for any non 
s-wave pairing state, the superconducting order on the 7 
sheet is weakly coupled to that on the q;,/3 sheets. ^"^^^^ 
This enables the physical picture of a multi gap super- 
conductor, where a full gap exists on the active band 
(either the 7 or the a, /3 bands) and the gap on the pas- 
sive band is driven by the active gap through the weak 
coupling between the two gaps. The gap on the passive 
band may or may not contain nodes depending upon the 
details of the coupling between the two bands. **'^^ This 
scenario is supported by the evolution of the specific heat 
with magnetic field where it is seen that a small magnetic 
field is sufficient to remove the gap in passive band.^ ' The 
identification of the active band has been a matter of de- 
bate. Specific heat measurements suggest that it is the 
7 band.^^ However, it has been recently argued that if 
the a, (3 bands are the active band, then it is possible 
that the edge currents are not large, *®'^^ which provides 
a reason for why they have not been observed (note that 
small edge currents may make it difficult to explain the 
fields seen by /zSR measurements'^"). Furthermore, un- 
like the 7 band, superconductivity in the a, /3 bands gives 
rise to an intrinsic Hall effect that naturally accounts for 
the observed polar Kerr effect. 

A striking feature which has been repeatedly observed 
in this material and is the subject of this paper, is vor- 
tex coalescence into clusters in various samples. ^'^"'^^'^^ 
The experimental works Refs. 9, 10, 22 interpreted this as 
originating from attractive intervortex interactions of un- 
known nature. Indeed the mechanisms for possible small 
intervortex attractions in single-component superconduc- 
tors cannot apply for this material. In single-component 
systems vortices have either purely repulsive (in type-2 
case) or purely attractive (in type-1 case) interaction at 
the level of Ginzburg-Landau theory. The situation can 
be more complex in microscopic models of weak-coupling 
type-2 superconductors with the Ginzburg-Landau pa- 
rameter K extremely close to the Bogomolnyi limit 1 / y/2. 
There, intervortex forces in Ginzburg-Landau theory are 
vanishingly small. The long-range intervortex forces are 
then determined by non-universal microscopic physics, 
which for certain materials can lead to tiny attractive 
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forces^^ which at long range can win over repulsive inter- 
action. However this mechanism can be ruled out in the 
case of Sr2Ru04 since the measurements of critical fields 
suggest that an estimate for a putative Ginzburg-Landau 
parameter is too far outside the regime where these ef- 
fects takes place: k « 2.6 (see e.g. Ref. 10). There ex- 
ists a different mechanism for intervortex attraction spe- 
cific for multicomponent systems. Multicomponent su- 
perconductors can have several coherence lengths ^a- In 
this case, the intervortex interaction can be long-range 
attractive and short-range repulsive because of the in- 
terplay of the multiple fundamental length scales of the 
theory. That is, a nonmonotonic interaction occurs if 
the London penetration length A falls between the coher- 
ence lengths ^1 < A < ^2-^^ "^^ This regime was recently 
termed "type-1.5" superconductivity^* since it features 
coexisting and competing type-I and type-II behaviors. 
The nonmonotonic intervortex forces originate there from 
the "double-core" structure of vortices due to multiple 
coherence lengths. The outer core (roughly speaking as- 
sociated to the length scale ^2) extends outside the flux 
carrying area. The overlap of the outer cores of vortices 
is responsible for the attractive intervortex forces. Yet in 
the type-1.5 regime the vortices have short range repul- 
sion (due to current-current and electromagnetic inter- 
action) and are thermodynamically stable. This should 
lead to vortex cluster formation in low magnetic field. 
Here we examine the role of the multiple superconducting 
degrees of freedom that exist in Sr2Ru04 . After it was 
suggested in Ref. 28 that type-1.5 regime is realized in 
MgB2 , a question was raised in the recent experimental 
work on Sr2Ru04^° whether or not the vortex coalescence 
in Sr2Ru04 originates via a type-1.5 scenario. Here we 
argue that for realistic choices of phenomenological pa- 
rameters, vortex clustering and type 1.5 behavior does 
occur in multiband chiral Ginzburg-Landau theories for 
Sr2Ru04 . 



We consider a two-band chiral Ginzburg-Landau the- 
ory which incorporates the multiband nature of the 
superconductivity. In accordance with the discussion 
above, the superconducting gap function is taken to be- 
long to the i?u representation of the tetragonal point 
group. This representation describes a chiral supercon- 
ductor. The four complex components of the order pa- 
rameter are related to the spin-triplet gap functions by 



da(k) = [i/'i 7x,a(k) -I- ?A^''Vy,a(k)]z whcrc the a la- 
bels one of the two bands and the functions fx,aO^) and 
/a,a(k) share the same symmetry properties as kx and ky 
under the symmetry operations of the point group -D4/1. 
In the units where ?i=l,c=l,m = l, the two-band 



situation can be modeled as follows: 
J" = |V X Ap + 

a=l,2 I 

+ 2Sa^aRe [{0x4"^)* Dy4"'' + {Dyil^[^^yDxi>i^^ 
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where ipj^ represent the superconducting components in 
the different bands, a = 1,2 denotes the band index, 
6 = 1,2 denotes the two different components of each 
condensate. Re[] stands for real part of the expression 
in brackets. Each component of a given condensate is 
a complex function '4' a = IV'al exp{j(/?a}- The gauge co- 
variant derivative is = V + ieA and B = V x A. 
The gauge coupling constant e is used to parametrize 
the penetration length of the magnetic field. This free 
energy is comprised of two free energies, one for each 
band, that are coupled through the vector potential A 
and by the parameter u. The free energies of each band 
are determined by weak coupling theory in the clean limit 
(both are reasonable assumptions for Sr2Ru04 ). No as- 
sumptions are made about the properties or geometry of 
the Fermi surface. The superconducting gap functions 
/j^ a(k) are assumed to be given by the Fermi velocity 
components Vxi}^) and Wy(k). This approximation pro- 
vides the correct momentum dependence of the gap for 
the limiting cases of a purely cylindrical Fermi surface 
and a perfectly square Fermi surface. The direct cou- 
pling of the bands through the sole parameter v is also 
justified within a weak coupling theory. The phenomeno- 
logical theory has nine parameters that require specifica- 
tion. Experiments place constraints on these parameters. 
In particular, one band is required to be passive and the 
other active. We assume that the 7-band is the active 
band and that the a/ j5 bands give rise to a passive gap. 
This assumption is not critical, the results will be similar 
if we assume the other possibility. This assump- 
tion implies ai < and a2 > 0. Since the realization of 
type-1.5 state will be very similar for ai > and a2 < 0, 
without loss of generality we report numerical investiga- 
tion of only the first scenario. 

Additionally, for magnetic fields applied along the c- 
axis, the free energy should reproduce the experimentally 
determined ratio of Hc2 and Hd. Finally, weak coupling 
calculations show that 71 > 1/3 and 72 < 1/3 and consis- 
tency with the small observed anisotropy of the in-plane 
upper critical field requires that both 71 and 72 are close 
to 1/3.33 
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Equations of motion of the gauge field defines the su- 
percurrent 



J- V J^"^ 



(2) 
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where the contribution of each component of a given con- 
densate is 
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where Im stands for imaginary part. 

Coalescence of vortex matter into clusters is inves- 
tigated numerically by minimizing the free energy (1) 
within a finite element framework provided by the 
FreeFemH — h library. '^^ We investigated vortex matter in 
the model (1) for parameters which give characteristic 
length scales which are close to experimental estimates 
for Sr2Ru04 . We report the case (ai,/3i) = (—10,10), 
7i — 0.35 and — 1. Parameters associated with the 
second condensate are (a2,/32) = (0.3, 1), 72 = 0.25 and 
82 = 1. The interband coupling is = 0.45 and the 
electric charge e = 1.8. We also investigated a range 
of similar parameters to ensure that they give a simi- 
lar picture (i.e. the parameter set does not correspond 
to any fine-tuned situation). For a general study how 
type-1.5 behavior is affected by interband coupling v see 
Refs. 25-27, and 29. Note that, the models like (1) can 
have a Skyrmionic phase, where the vortices are unstable 
against a decay to Skyrmions."^^ We obtained stable vor- 
tex clusters in various type-1.5 regimes which are stable 
against a decay into Skyrmions. Fig. 1 shows a numeri- 
cal solution for one such vortex cluster. We also find that 
Skyrmions can still exist in those regimes as metastable 
topological excitations, which are more energetically ex- 
pensive than vortices. 

Note that the experiments^'^"'^^ were performed at 
T <C Tc- A microscopic approach is therefore required to 
get all aspects of the physics quantitatively correct since 
the Ginzburg-Landau model is quantitatively correct 
only for T near to T^. However, under certain circum- 
stances, a phenomenological Ginzburg-Landau model 
can qualitatively and quantitatively describe the long- 
range intervortex forces in the type-1.5 regime even at rel- 
atively low temperatures. This was shown in a two-band 
system with not too strong interband coupling and one 
passive band.'^'^ The physics that a GL-based approach 
fails to describe in the low temperature regime, is pri- 
marily associated with shorter-length scales such as the 
counterpart of Kramer-Pesch effect. Since in our study 
we are only interested in long-range intervortex forces 
(i.e. we do not consider the high-field regimes, where in- 
tervortex distance is short), the phenomenological GL 
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Figure 1. (Color online) - A vortex cluster obtained nu- 
merically by energy minimization from a dilute system of 
seven vortices. Displayed quantities are the magnetic field 
on (A) and the total supercurrent (2) on (B). The den- 
sities of the first {resp. second) component of the conden- 
sate 'ip''^\ namely jV'i^'P {resp. \ip2^^\^), are shown on (C) 
{resp. (D)). The panels (E) {resp. (F)) show the densities of 
the first {resp. second) component of the condensate i/j'^-' , that 

is 1^1^' 1^ {resp. 1V'2^'P)- The figure demonstrates existence of 
several length scales associated with the density variations. 
The panels (E) {resp. (F)) show extended cores, the overlap 
of these cores leads to attractive intervortex forces. 



approach still provides a qualitatively correct picture of 
that physics even at relatively low temperatures. 

The experiments in Refs. 22, 10 and 23 obtain rather 
similar intervortex distances within the clusters. The 
typical intervortex distance is larger than their estimates 
of the London penetration length. Intervortex attraction 
is indeed achievable in type-1.5 regime at such length 
scales. The minimum of the intervortex interaction po- 
tential is determined not only by the length scales A, ^1, ^2 
but also by non-linear effects. Thus in a type-1.5 regime 
intervortex distance in a cluster could be substantially 
larger than the London penetration length scale. For 
example, the intervortex potential shown on Figure 2 
in Ref. 26 corresponds to such a situation. It should 
be noted that all experiments observing vortex clusters 
(Refs. 9, 10, 22), are scanning SQUID or Hall experiments 
which probe magnetic field at a small distances over the 
sample's surface. It should be kept in mind that such a 
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surface probe could in general overestimate the position 
of the minimum of the interaction potential, for vortices 
with nonmonotonic interactions. Indeed, near the sur- 
face, the long-range intervortex forces can be altered by 
the electromagnetic repulsion caused by the stray fields 
outside the sample (schematically shown in Fig. 2). Sim- 
ilarly, stray field physics affects the structure of normal 
domains in ordinary type-1 superconductors.'^*' 

stray magnetic field 




Figure 2. Even if vortices have attractive forces in the bulk, 
their segments near the surface can have long-range contribu- 
tions to repulsive interaction due to the effects of demagneti- 
zation fields outside the sample. This is schematically shown 
on this figure. Thus a surface probe of a bulk sample can 
overestimate the position of the minimum in the intervortex 
potential. Also in case of a thin film, the stray fields can make 
intervortex distance larger. 

The simplest origin of inhomogeneous vortex distribu- 
tions are the pinning effects. And in fact in some of the 
samples, existence of preferential pinning areas were iden- 
tified in the work of Curran et. al.^^^'^'^ see also remark"^*. 
However, in the experimental papers^^ it was argued that 
the observed vortex clustering is not related to pinning. 
The layered structure of Sr2Ru04 leads to vortex stripes 
formation when in-plane field is applied. Therefore it is 
possible to move vortices by applying the in-plane com- 
ponent of the external magnetic field. Thus it is possi- 
ble to assess this way, if the vortex clusters are artifacts 
of some local pinning landscape. In the Ref. 22 it was 
argued that the mobility of the vortex clusters in their 
samples is inconsistent with clustering due to a pinning 
scenario. 

The distribution of the intervortex distance in several 



samples with a relatively small number of vortices was 
analyzed in.^'^ They found that it did not exhibit a clear 
peak at certain preferred intervortex distance, which was 
interpreted as being inconsistent with the type- 1.5 sce- 
nario. We argue however that it does not necessarily 
contradict the type- 1.5 scenario, for the following rea- 
sons. The shape of the vortex clusters in a certain sub- 
set of type-1. 5 regimes is rather substantially affected by 
non-pairwise contributions to intervortex forces^° (and 
also by dynamic and entropic aspects associate with it). 
This can affect the distribution of intervortex distances. 
Since the multiband model of Sr2Ru04 has four com- 
ponents, the non-pairwise contributions to intervortex 
forces could be relatively significant and produce vari- 
ation in the shapes of vortex clusters. Due to these ef- 
fects, it can be more difficult to get a distinguishable 
double-peak structure from vortex distributions in small 
samples like those studied in Ref. 23. Large samples with 
a larger number of vortices can yield a more conclusive 
answer. The absence of collapse to a single vortex clus- 
ter noticed in Refs. 22 and 10, also naturally arises in 
type- 1.5 regime. Again it can originate in dynamic and 
entropic reasons and be enhanced by non-pairwise con- 
tributions to intervortex forces. It can also be a result 
of demagnetization (stray fields) effects. 

In conclusion, we have shown that vortex coalescence 
occurs in a realistic phenomenological model describ- 
ing the multicomponent superconductivity in Sr2Ru04 . 
This provides an explanation for the vortex coalescence 
that was reported by several experimental groups. These 
experiments can be interpreted as demonstrating that 
Sr2Ru04 is a type-1. 5 superconductor: i.e. it possess 
several coherence length: some longer than the London 
penetration length and some shorter. In order to firmly 
establish this interpretation, further experimental studies 
on the distribution of vortices are required. 

We thank Peter Curran, Clifford Hicks, Simon Bend- 
ing, Kam Moler and Victor Moshchalkov for useful cor- 
respondence. This work was supported by Knut and Al- 
ice Wallenberg Foundation through the Royal Swedish 
Academy of Sciences, Swedish Research Council and by 
the NSF CAREER Award No. DMR-0955902. DFA 
is supported by NSF grant DMR-0906655. The com- 
putations were performed on resources provided by the 
Swedish National Infrastructure for Computing (SNIC) 
at National Supercomputer Center at Linkoping, Swe- 
den. 



^ A. P. Mackenzie and Y. Maeno, Rev. Mod. Phys. 75, 657 
(May 2003) 

^ K. Ishida, H. Mukuda, Y. Kitaoka, K. Asayama, Z. Q. 
Mao, Y. Mori, and Y. Maeno, Nature 396, 658 (Dec. 1998) 

^ H. Murakawa, K. Ishida, K. Kitagawa, Z. Q. Mao, and 
Y. Maeno, Phys. Rev. Lett. 93, 167004 (Oct. 2004) 
G. M. Luke, Y. Fudamoto, K. M. Kojima, M. I. Larkin, 
J. Merrin, B. Nachumi, Y. J. Uemura, Y. Maeno, Z. Q. 



Mao, Y. Mori, H. Nakamura, and M. Sigrist, Nature 394, 

558 (Aug. 1998) 
^ J. Xia, Y. Maeno, P. T. Beyersdorf, M. M. Fejer, and 

A. Kapitulnik, Phys. Rev. Lett. 97, 167002 (Oct. 2006) 
^ K. D. Nelson, Z. Q. Mao, Y. Maeno, and Y. Liu, Science 

306, 1151 (2004) 
^ F. Kidwingira, J. D. Strand, D. J. Van Harlingen, and 

Y. Maeno, Science 314, 1267 (2006) 



5 



® M. Matsumoto and M. Sigrist, J. Phys. Soc. Jpn. 68, 994 
(1999) 

^ P. G. Bjornsson, Y. Maeno, M. E. Huber, and K. A. Moler, 

Phys. Rev. B 72, 012504 (Jul. 2005) 
^° C. W. Hicks, J. R. Kirtley, T. M. Lippman, N. C. Koshnick, 
M. E. Huber, Y. Maeno, W. M. Yuhasz, M. B. Maple, and 
K. A. Moler, Phys. Rev. B 81, 214501 (Jun. 2010) 
J. Jang, D. G. Ferguson, V. Vakaryuk, R. Budakian, S. B. 
Chung, P. M. Goldbart, and Y. Maeno, Science 331, 186 
(Jan. 2011) 

D. F. Agterberg, Phys. Rev. Lett. 80, 5184 (Jun. 1998) 
C. Bergemann, A. P. Mackenzie, S. R. Julian, D. Forsythe, 
and E. Ohmichi, Advances in Physics 52, 639 (2003) 
" D. F. Agterberg, T. M. Rice, and M. Sigrist, Phys. Rev. 
Lett. 78, 3374 (Apr. 1997) 

M. E. Zhitomirsky and T. M. Rice, Phys. Rev. Lett. 87, 
057001 (Jul. 2001) 

J. F. Annett, G. Litak, B. L. Gyorffy, and K. I. Wysokihski, 
Phys. Rev. B 66, 134514 (Oct. 2002) 

K. Deguchi, Z. Q. Mao, H. Yaguchi, and Y. Maeno, Phys. 
Rev. Lett. 92, 047002 (Jan. 2004) 

S. Raghu, A. Kapitulnik, and S. A. Kivelson, Phys. Rev. 
Lett. 105, 136401 (Sep. 2010) 

Y. Imai, K. Wakabayashi, and M. Sigrist, Phys. Rev. B 85, 
174532 (May 2012) 
^" P. E. C. Ashby and C. Kalhn, Phys. Rev. B 79, 224509 
(Jun. 2009) 

2^ E. Taylor and C. Kallin, Phys. Rev. Lett. 108, 157001 
(Apr. 2012) 

^'^ V. O. Dolocan, C. Veauvy, F. Servant, P. Lejay, K. Has- 
selbach, Y. Liu, and D. Mailly, Phys. Rev. Lett. 95, 
097004 (Aug. 2005)K. Hasselbach, V. Dolocan, P. Lejay, 
and D. Mailly, Physica C 460-462, Part 1, 277 (2007) 

2^ P. J. Curran, V. V. Khotkevych, S. J. Bending, A. S. Gibbs, 
S. L. Lee, and A. P. MacKenzie, Phys. Rev. B 84, 104507 
(Sep. 2011) 

2" A. Jacobs, J. Low Temp. Phys. 10, 137 (1973)M. C. Leung 
and A. E. Jacobs, ibtd. 11, 395 (May 1973)M. C. Leung, 
ibid. 12, 215 (Jul. 1973)U. Klein, ibid. 69, 1 (1987) 

'^^ E. Babaev and J. M. Speight, Phys. Rev. B 72, 180502 
(2005)E. Babaev, J. Carlstrom, and M. Speight, Phys. Rev. 
Lett. 105, 067003 (Aug. 2010) 



^"^ J. Carlstrom, E. Babaev, and M. Speight, Phys. Rev. B 
83, 174509 (May 2011) 

M. Silaev and E. Babaev, Phys. Rev. B 84, 094515 (Sep. 
2011) 

V. Moshchalkov, M. Menghini, T. Nishio, Q. H. Chen, 
A. V. Silhanek, V. H. Dao, L. F. Chibotaru, N. D. Zhi- 
gadlo, and J. Karpinski, Phys. Rev. Lett. 102, 117001 
(Mar. 2009)T. Nishio, V. H. Dao, Q. Chen, L. F. Chib- 
otaru, K. Kadowaki, and V. V. Moshchalkov, Phys. Rev. 
B 81, 020506 (Jan. 2010) 

M. Silaev and E. Babaev, Phys. Rev. B 85, 134514 (Apr. 
2012) 

J. Carlstrom, J. Garaud, and E. Babaev, Phys. Rev. B 84, 
134515 (Oct. 2011) 

V. H. Dao, L. F. Chibotaru, T. Nishio, and V. V. 
Moshchalkov, Phys. Rev. B 83, 020503 (Jan. 2011) 

3^ R. Geurts, M. V. Milosevic, and F. M. Peelers, Phys. Rev. 
B 81, 214514 (Jun. 2010) 

^3 D. F. Agterberg, Phys. Rev. B 64, 052502 (Aug. 2001) 
F. Hecht, O. Pironneau, A. Le Hyaric, and K. Ohtsuka, 
The Freefem++ manual (2007) www.freefem.org 
The integration domain, is discretized using an algorithm 
based on Delaunay-Voronoi triangulation and each func- 
tion is decomposed on a piecewise quadratic basis. The 
discrete version of Eq. (1) is then minimized with respect 
to all degrees of freedom, for a starting configuration with 
an initial guess for A'^ vortices using a Non-Linear Conju- 
gate Gradient scheme. 

^5 J. Garaud and E. Babaev (Jan. 2012), arXiv:1201.2946 

^® P. de Gennes, Superconductivity of metals and alloys 
(Perseus Books, 1999) 

P. J. Curran and S. J. Bending, Private communication 
Note that in chiral superconductors vortex pinning can 
arise not only from sample inhomogeneities but also due 
to vortex interaction with pinned domain walls between 
different chiral domains'*^'''^ . 

M. Ichioka, Y. Matsunaga, and K. Machida, Phys. Rev. B 
71, 172510 (May 2005)V. Vakaryuk, ibid. 84, 214524 (Dec. 
2011)D. G. Ferguson and P. M. Goldbart, Phys. Rev. B 84, 
014523 (Jul. 2011) 



